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The Link Between the Formation Rates of Clusters and Stars in 

Galaxies 
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ABSTRACT 

The goal of this paper is to test whether the formation rate of star clnsters 
is proportional to the star formation rate (SFR) in galaxies. As a hrst step, we 
present the mass fnnctions of compact clusters younger than 10 Myr in seven 
star-forming galaxies of diverse masses, sizes, and morphologies: the Large and 
Small Magellanic Clouds, NGC 4214, NGC 4449, M83, M51, and the Anten¬ 
nae. These cluster mass functions (CMFs) are well represented by power laws, 
dN/dM oc M^, with similar exponents fd = —1.92 ± 0.27, but with amplitudes 
that differ by factors up to ~10^, corresponding to vast differences in the sizes of 
the cluster populations in these galaxies. We then normalize these CMFs by the 
SFRs in the galaxies, derived from dust-corrected Ha luminosities, and hnd that 
the spread in the amplitudes collapses, with a remaining rms deviation of only 
cr(logA) = 0.2. This is close to the expected dispersion from random uncertain¬ 
ties in the CMFs and SFRs. Thus, the data presented here are consistent with 
exact proportionality between the formation rates of stars and clusters. However, 
the data also permit weak deviations from proportionality, at the factor of two 
level, within the statistical uncertainties. We hnd the same spread in amplitudes 
when we normalize the mass functions of much older clusters, with ages in the 
range 100 to 400 Myr, by the current SFR. This is another indication of the 
general similarity among the cluster populations of different galaxies. 

Subject headings: galaxies: individual (LMC, SMC, NGC 4214, NGC 4449, M83, 
M51, Antennae) — galaxies: star clusters — stars: formation 


1. Introduction 

Stars form together in clusters and associations, which in turn form in the densest 
parts—the clumps—of molecular clouds (Lada & Lada 2003; McKee & Ostriker 2007). There 
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is some debate about the exact definition of a cluster, e.g., whether this term refers to all 
dense stellar aggregates or only those that are gravitationally bound or free of gas. Nev¬ 
ertheless, there appears to be unanimous agreement that stars mostly form in aggregates 
of some sort and rarely in isolation. This general picture implies that the formation rates 
of stars and clusters should be proportional to each other. The pnrpose of this paper is to 
present a direct test of this proportionality. 

Our approach is to compare the mass functions, ■^(M) = dN/dM, of recently formed 
star clusters (with ages r < 10 Myr) in different galaxies before and after normalizing by 
the star formation rates in the galaxies. This is valid because the mass functions of yonng 
clusters in different galaxies have similar power-law shapes, oc with {3 ~ —2, over 

large intervals of mass, roughly 10^ Mq ^ M ^ 10® Mq (Fall & Chandar 2012 and references 
therein). The mass fnnctions are unknown at lower masses, as a result of selection limits 
in the available samples, and are nncertain at higher masses as a result of small number 
statistics. This is the reason that we compare the mass functions, rather than the total 
nnmbers or total masses of yonng clusters (the integrals of ip{M) and M'lpiyM) respectively 
over all M), with the star formation rates. 

In this work, we use the term “cluster” for any concentrated aggregate of stars with a 
density mnch higher than that of the snrronnding stellar field, whether or not it is gravita¬ 
tionally bound. It is not possible to tell from images alone which clusters are gravitationally 
bonnd (have negative energy) and which are nnbound (have positive energy), and N-body 
simnlations show that unbound clusters retain the appearance of bound clusters for many 
(~10) crossing times (Banmgardt & Kronpa 2007). While spectra might help in principle, 
they are seldom available and are often contaminated by non-virial motions (stellar winds, 
binary stars, etc). 

In the next section (Section 2), we present the catalogs and mass fnnctions in seven 
galaxies with well-stndied popnlations of clusters. In the following section (Section 3), we 
estimate the star formation rates, determine the CMF/SFR relations, present a method for 
comparing this relation between galaxies, and compare the residuals with various galaxy 
properties. In the last section (Section 4), we snmmarize some implications of onr results 
for the formation and disrnption of the clusters. 


2. Cluster Catalogs and Mass Functions 


Our study focuses on seven nearby galaxies: LMC, SMC, NGC 4214, NGC 4449, M83, 
M51, and the Antennae. In this section, we collect and describe the star cluster catalogs. 



and summarize the method used to estimate the masses and ages of the clusters. We then 
present the mass functions of the clusters in each galaxy. 


2.1. Star Cluster Samples 

We include nearby galaxies with published star cluster catalogs that cover at least 
50% of the optically luminous portions of the galaxy and have photometry in at least four 
optical passbands including the U band (which is critical for disentangling the effects of age 
and extinction in the colors of young clusters). Our sample includes galaxies of different 
types: irregular, spiral, and merger, which span a large range in distance, luminosity, color, 
inclination, and star formation rate. This sample, while small, is reasonably representative 
of nearby star-forming galaxies in general. The observations, data, and cluster selection 
criteria for each galaxy are summarized below. 

Large and Small Magellanic Clouds: We use the catalogs published by Hunter et ah 
(2003), where clusters were selected by visual examination of candidates (from previously 
published catalogs). The clusters were distinguished from the surrounding distribution of 
held stars, resolved with respect to a single star, but do not include low density H ll regions. 
The observations cover approximately 70% and 90% of the recent star formation in each 
galaxy, as traced by Ha emission. This resulted in a total of 854 (239) clusters in the LMC 
(SMC) with published UBVR photometry. 

NGC 4214: This dwarf irregular galaxy, at a distance of 3.1 Mpc (Dopita et al. 2010), 
was observed in a single pointing with the WFC3 on-board the Hubble Space Telescope 
{HST) as part of program 11360 (PI: O’Connell). Clusters were automatically selected to 
be compact sources broader than the PSF followed by elimination of close pairs of stars, 
which resulted in a total of 334 candidate clusters. Aperture photometry was performed on 
UBVIHa images in a manner similar to that described in Chandar et al. (2010) for sources 
in the nearby spiral galaxy M83. 

NGC 4449: Rangelov et al. (2011) published a catalog of 129 clusters in this dwarf 
starburst galaxy (located at a distance of 3.82 ± 0.27 Mpc; Annibali et al. 2008) based 
on HST observations. They selected cluster candidates automatically using size criteria, 
followed by a visual inspection to eliminate contaminants such as background galaxies and 
chance superpositions. Two helds within NGC 4449 were observed with the ACS/WFC in 
the BVIHa bands, and two in the U hlter with the WFPC2 camera, giving near-complete 
coverage of the optically luminous portions of the galaxy. 
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M83: Two fields within M83, a grand-design spiral galaxy located at a distance of 
4.5 Mpc (Thim et al. 2003), were observed as part of program GO-11360. Several cluster 
catalogs were produced for each held, according to the selection methodology described in 
Chandar et al. (2014). Here, we adopt the manual cluster catalogs from that work. We 
also use data from hve additional helds in M83 that were observed with the HST /WFC3 as 
part of program GO-12513, bringing the coverage of the optically luminous portion of M83 
to ~60%. New cluster catalogs have been produced for these helds, and will be presented in 
Whitmore et al. (in prep). Our hnal catalog contains 3186 star clusters. 

M51: Most of the optically luminous portions of this grand-design spiral galaxy, located 
at a distance of 8.2 Mpc, were observed in six helds taken with the HST/ACS camera 
in BVIHa, and in 8 overlapping helds with the WFPG2 camera in the U band. Just 
as for NGG 4449, star clusters were selected using automated criteria, followed by visual 
examination. This resulted in a catalog of 3812 clusters (Ghandar et ah, in prep). 

Antennae: The merging Antennae galaxies, located at a distance of 21 Mpc (Schweizer 
et al. 2008), were observed with the AGS/WFG on the HST as part of program GO-10188, 
in the BVIHa hlters (Whitmore et al. 2010). Additional observations were taken in the 
U band with the WFG3 camera as part of program GO-11577. Glusters were selected to be 
objects brighter than My = —9, a restriction that eliminates practically all individual stars. 

The largest source of uncertainty in these catalogs is in the selection of clusters with 
ages r < 10 Myr. For the Magellanic Glouds, it is possible that clusters that are both 
very young and/or have low stellar density were excluded by the selection criteria. For the 
more distant galaxies, very young clusters might be excluded because they are very compact 
and their prohles are indistinguishable from individual stars, but on the other hand, chance 
superpositions of unassociated young stars might sometimes be selected as young clusters. 
The selection of compact clusters older than 10 Myr is more robust (as discussed in Ghandar 
et al. 2014), because by this age clusters have typically moved away from their crowded birth 
sites. 


2.2. Mass and Age determinations 

For each cluster in each galaxy, the measured magnitudes, in the hlters listed above, 
were compared with predictions from the Bruzual & Gharlot (2003, 2010) stellar population 
models, in order to estimate their age, extinction, and mass. 

We hrst estimate the age r and the extinction Ay by performing a fif comparing 
observed and predicted magnitudes, assuming a Ghabrier (2003) IMF, and a Galactic-type 
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extinction law (Fitzpatrick 1999). The best-fit values of r and Ay are those that minimize: 

= • ( 1 ) 

A 

where and are the observed and model magnitudes respectively, and the sum 

runs over all hlters listed above. The weight factors in the formula for are taken to be 
W\ = [a\ + (0.05)^]“^, where a\ is the photometric uncertainty. 

The mass of each cluster is estimated from the extinction-corrected Id-band magnitude, 
the distance listed in Table 1, and the mass-to-light ratio {M/Ly) predicted by the models at 
the htted ages r. The models assume that the stellar IMF for each cluster is fully sampled. 
While the colors of clusters with masses below Mq begin to spread around the mean 
values because the upper end of the stellar IMF is not fully populated, the resulting mass 
distributions, down to at least x 10^ Mq, are very similar regardless of whether fully or 
stochastically sampled models are used (Fouesneau et ah 2012). 

The largest uncertainties in the mass estimates comes from uncertainties in the age 
estimates, which are typically on the order of 0.3 in logr, corresponding to a factor of ~2 
in r. These translate to approximately similar uncertainties of 0.3 and Ri2 in logM and 
M. There are also two sources of systematic uncertainty. Any error in a galactic distance 
translates directly into a fractional change in the masses of all clusters within that galaxy. 
The derived masses of the clusters, but not their ages, also depend on the IMF assumed in 
the stellar population models. For example, if we had assumed a Salpeter (1955) instead of 
Chabrier IMF, the estimated mass of each cluster would increase by a nearly constant ~40%. 
Similarly, we hud that a mix of stellar IMFs between Chabrier and Salpeter would increase 
the estimated masses by a nearly constant fractional amount between 0 and 40% (depending 
on the exact mix of IMFs). Such fractional changes in M, corresponding to shifts in logM, 
would not alter the shape of the cluster mass function (in log-log plots). 


2.3. Cluster Mass Functions 

In the left-hand panels of Figures 1 and 2, we show the mass functions determined for the 
clusters in all seven of our sample galaxies in two different intervals of age: recently formed 
(r < 10 Myr) clusters, which are the main focus of this paper, and also intermediate-age 
(r = 100-400 Myr) clusters as an interesting check. An equal number of clusters have been 
placed in each bin, from 3 to 7, depending on the available sample size (resulting in variable 
bin widths). The low end of the mass function in each galaxy is set by the luminosity limit 
of each catalog, while the high end of these distributions is not constrained. The amplitudes 
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of the cluster mass functions directly reflect differences in the sizes of the cluster populations 
among the galaxies. There is clearly a large range, approximately a factor of 10^, in these 
amplitndes. 

The merging Antennae system has the highest nnmber of stellar clnsters in onr sample 
(in both age intervals). It is followed by the spirals M51 and M83, then the Inminons star¬ 
forming irregnlar galaxies NGC 4449 and the LMC, with the low Inminosity irregnlar galaxies 
NGC 4214 and the SMG having the fewest detected clusters. The high end of the clnster 
mass fnnctions also generally correlate with the total nnmber of clusters, in the sense that 
the galaxies with the most clusters tend to have the most massive clusters. 

The observed mass functions are well represented by featureless power laws: dN/dM oc 
M^. Thus, we perform linear hts to log((iA^/(ilogM) = {/3 + l)logM-|- const. The best- 
£t exponents and their formal la errors are listed in Table 1. Despite the very different 
types of galaxies inclnded in onr sample, different cluster selection criteria, filter sets, etc., 
the resulting exponents have a fairly narrow range, /d = —1.92 ± 0.27, indicating that the 
mass fnnctions of yonng star clusters in nearby star-forming galaxies have an approximately 
“nniversal” shape. 

Uncertainties in the GMF are dominated by the selection of the clusters. For clusters 
younger than 10 Myr, we estimate that the nncertainties in the amplitude of the GMF could 
be as large as ^0.1 based on comparisons of the cluster catalogs from different gronps, which 
differ at the 30-40% level in this age range (e.g., Bastian et ah 2012; Ghandar et al. 2014). 
The nncertainty for 100-400 Myr clnsters is smaller, ronghly Ri0.04 because cluster catalogs 
in this age range only differ by ~10% between gronps. 


3. Star Formation Rates and Normalized Mass Functions 
3.1. Star Formation Rates 

We determine the star formation rate for each galaxy in our sample from its Ha line 
emission. Most of this emission (~90%) is contribnted by stars yonnger than 10 Myr (Ken- 
nicntt & Evans 2014). Therefore, Ha emission measures the SFR over the same age range 
as the clusters in onr yonnger snbsample. We use the most recent Ha flnx measnrements 
listed in NED and corrected for contamination by [N ll] emission, and assnme the distances 
given in colnmn 3 of Table 1 when converting to Ha Inminosities. 

The Ha emission is attennated by dnst. To correct for this attennation, we use the 
measured 25 /rm flux (available for all galaxies in onr sample), and the recipe given in 
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Table 2 in Kennicutt & Evans (2014): 

L[Ha)coTT = + 0.020L(25/i?7i) . (2) 

We then adopt their Equation (12) and Table 1 to determine the star formation rate: 

log SFR(M 0 /yr) = log L(Ha)corr - 41.27 . (3) 

This calibration is based on the STARBURST99 models (Leitherer et al. 1999) assuming solar 
metallicity, a Kroupa IMF (similar to the Chabrier IMF in the Bruzual & Chariot models), 
and a constant rate of star formation. The resulting SFRs are compiled in column 5 of 
Table 1. 

There are several inherent sources of uncertainty in Ha-based star formation rate deter¬ 
minations, including flux measurements, corrections made for attenuation by dust, variations 
in the star formation rate, uncertainties in the calibration used to convert the measured 
fluxes to SFR as discussed in the review by Kennicutt & Evans (2014), and the leakage 
of Lyman continuum photons from the parent galaxy. Estimates of these uncertainties are 
listed in Table 2. For our sample specihcally, we hud that the published Ha and 25 yum flux 
measurements (compiled in NED) used to determine the SFRs, typically differ by ~20%, 
although one galaxy has published values that differ by ~50%. Different methods to correct 
for dust extinction introduce uncertainties of ~20-30%, based on a comparison between a 
linear combination of the Ha and 24 /rm luminosities with Balmer-decrement corrected Ha 
luminosities (Figure 3, Kennicutt & Evans 2014), although the uncertainties can be higher 
for dusty galaxies like the Antennae. We estimate that the uncertainty in the fractional 
coverage of each galaxy is ~10%. Lyman continuum photons can escape from their parent 
galaxy, further increasing the uncertainty in the SFR determinations, although the escape 
fraction can vary strongly from one galaxy to another and is therefore difficult to quantify. 
Taken together, the random uncertainties related to the SFR determinations and coverage 
correspond to at least RiO.2 in the logarithm of the CMF/SFR relation, not including uncer¬ 
tainties related to the cluster mass functions. Systematic uncertainties may also be present, 
since stellar population modeling shows that the calibrations likely overestimate the SFR for 
lower metallicity populations, because the ionizing luminosity increases by ~0.4 ±0.1 dex 
for a tenfold decrease in the metallicity (e.g., Kennicutt & Evans 2014; Smith, Norris, & 
Crowther 2002; Raiter, Schaerer & Fosbury 2010). 

Our Ha-based star formation rates agree reasonably well with previous determinations 
published in the literature (compiled in Column 6 of Table 1). Lee et al. (2009) determined 
Ha-based SFRs for NGC 4214, NGC 4449, M83, and M51, but used a different method to 
correct for dust and a different calibration to convert Ha luminosity to SFR (taken from 
Kennicutt 1998). James et al. (2008) estimated Ha-based SFRs for the Large and Small 



Magellanic Cloud using the prescription given by Kennicutt (1998), and Zhang et ah (2001) 
measured the Ha line intensity from HST /WFPC2 images and applied the prescription from 
Kennicutt (1998) to determine the SFR of the Antennae. In Column 7 of Table 1, we list 
previously published star formation rates determined for each galaxy using various non-Ha- 
based methods. Overall, we hnd that our SFR determinations agree with published values 
to somewhat better than a factor of two, which is consistent with the uncertainties found by 
Lee et ah (2009) by comparing far ultraviolet and Ha-based SFRs for a much larger sample 
of galaxies. 


3.2. The CMF/SFR Relation 

The right panel of Figure 1 presents the main result of this paper. The same cluster 
mass functions as in the left panel are shown, but these distributions have now been divided 
by the star formation rate (listed in Column 5 of Table 1) of the host galaxy within the 
area covered by each cluster catalog: CMF/SFR. The mass functions of star clusters in 
these very different galaxies all have very similar amplitudes when normalized by the star 
formation rate of their host galaxies. 

In Figure 2, we again show the observed and SFR-normalized cluster mass functions, 
but now for clusters with ages between 100 and 400 Myr. Evidently, the mass functions 
of these intermediate-age clusters also lie very close to one another in the vertical direction 
after they have been normalized by the current star formation rate of the host galaxy. 
This is remarkable because the SFR-normalized mass function of intermediate-age clusters 
can be affected by variations in the past formation and disruption rates of the clusters. The 
similarity of these CMFs indicates that the SFRs were relatively steady for the past 400 Myr, 
and that the cluster disruption rates were similar among the galaxies in our samples. 

In order to quantify the observed scatter, we ht the CMF/SFR for each galaxy by the 
function: 

dN/dM = Ax SFR x (M/10^ , (4) 

where we have hxed the power-law exponent to —1.9 (the mean of all ft values) for all galaxies. 
We use our Ha-based SFRs from Table 1, adjust them for partial coverage (if needed), then 
normalize the distributions at 10^ Mq, which is about the middle of the cluster log mass 
range observed across the galaxies. The coefficient A measures the proportionality between 
the cluster and star formation rates. The dispersion a (log A) in the best ht values of log A, 
quantihes the scatter in the amplitudes and hence in the CMF/SFR relation among galaxies. 

We demonstrate our technique in Figure 3. The SFR-normalized cluster mass functions 
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are shown in two intervals of age: r < 10 Myr (left panel) and r = 100-400 Myr (right 
panel). The lines show the best £t from equation (4), and the intercept values A are the 
values where each fit intersects the dotted line. For the seven galaxies studied here, we 
find that (T(log A) = 0.21 for clusters with ages r < 10 Myr. The scatter for intermediate- 
age (100-400 Myr) clusters is similar, with (j(log A) = 0.20. This scatter is close to the 
dispersion expected from random uncertainties in the CMFs and the SFRs, as discussed in 
Sections 2.3 and 3.1. 

The data presented here are consistent with exact proportionality between the formation 
rates of stars and clusters. Nevertheless, within the uncertainties, the data also allow for 
weak deviations, at the approximately factor of two level. In Figure 4, we search for trends in 
the residuals of A with host galaxy properties, including SFR, luminosity, color, inclination, 
maximum rotational velocity of the gas, and the exponent of the cluster mass function /3. 
The nonparametric Spearman correlation coefficients, rs and the corresponding probabilities 
p(> f's) associated with each plot are compiled in Table 3. Values for p(> rs) should be 
less than 0.05 to have greater than 95% confidence that a correlation is real. For the young 
(r < 10 Myr) clusters, we do not observe any statistically significant correlations. For the 
intermediate-age (100-400 Myr old) clusters, there are weak positive trends between the 
residuals of log A and galactic SFR, luminosity, and B—V color, but only the one with B—V 
is statistically significant (at the 95% level). Possibly, these trends may reflect differences 
in the metallicities of the galaxies. In any case, a larger sample of galaxies is needed to 
determine whether or not such trends are real (e.g., the LEGUS project; Calzetti et al. 
2015). We find no statistically significant correlation for clusters in either age range between 
the residuals of log A and the SFR per unit area or SFR density (not shown), which is often 
used to compare cluster populations among galaxies. 


3.3. Comparison with Previous Work 

The CMF/SFR statistic studied here is closely related to the cluster formation efficiency 
F, dehned as the fraction of stellar mass formed in bound clusters. Estimates of F in 
numerous studies have been used to draw important conclusions about the properties of 
star-forming molecular clouds in different galaxies and environments. Since, as discussed 
in the Introduction, it is not possible to distinguish bound from unbound clusters based on 
images alone, F should be directly proportional to our CMF/SFR statistic. 

F has been estimated for a number of galaxies, using various methods, and large differ¬ 
ences are claimed from one galaxy to another. All of these studies integrate over the CMF in 
order to determine the total mass in clusters, but use different assumptions, extrapolations. 
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age and mass ranges, and methods to accomplish this (e.g., Goddard et ah 2010; Adamo 
et ah 2011; Silva-Villa & Larsen 2011; Silva-Villa et ah 2013). The compiled results (e.g., 
Kruijssen 2012) show a strong increase in T, by a factor of ~15, for galaxies with increasing 
star formation rate densities, from Ril0“^ to ~ a few Mq yr“^ kpc“^ (corresponding to the 
~3 dex range of SFRs from ~0.05 to ~50 Mq yr“^). The corresponding relation between T 
and SFR (which has a similar slope when transformed from the original SFR density given 
in Kruijssen 2012) is shown as the dashed line in the top-left panel of Figure 4. This trend 
is much stronger than our CMF/SFR relation for young (r < 10 Myr) and intermediate-age 
(100-400 Myr) clusters. 

This may be due to the fact that estimates of F are sensitive to the specihc assumptions 
and extrapolations that are madeji] and are thus prone to large errors. For example, God¬ 
dard et ah (2010) found a signihcantly higher value of F for clusters in the Antennae when 
compared with those in the Magellanic Glouds, while we find similar GMF/SFR relations 
for these three galaxies. Our results are based on a larger cluster catalog for the Antennae, 
but the same ones as used in Goddard et al. (2010) for the LMG and SMG. The dotted line 
in the top-left panel of Figure 4 shows that the theoretical prediction for F versus SFR from 
Kruijssen (2012) is much steeper than the observations, and therefore disagrees even more 
strongly with our results for the GMF/SFR statistic. 

We have presented a robust new method for comparing the fraction of stars that form 
in clusters, one that only relies on comparing the amplitudes of the GMF normalized by the 
SFR in galaxies. The parameter F has been extensively used in the literature, but requires 
additional assumptions and extrapolations that degrade its reliability. We also advocate 
that, as done here, the SFR rather than the more uncertain SFR density be used when 
comparing results among galaxies. There is no standard method for estimating the star¬ 
forming area within galaxies, and we find that the uncertainties in area alone can be as large 
as a factor of a few (see also discussion in Adamo et al. 2011). 


4. Conclusions and Implications 

The main conclusion of this paper is that star and cluster formation rates are propor¬ 
tional to each other in a diverse but representative sample of nearby galaxies. There are 
only minor deviations from exact proportionality in our sample (with an rms scatter less 


^In particular, in this approach the extrapolated number of clusters depends sensitively on the ht and 
uncertainties in the power-law index /3 and also on the exact mass and age ranges that are adopted. Small 
uncertainties in any of these values can drastically change the results. 
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than a factor of 2) in two different age ranges, and no trends are observed in the residnals 
of CMF/SFR with galaxy property for young (r < 10 Myr) clusters. For the intermediate- 
age (100-400 Myr) clusters, weak trends with the star formation rate and luminosity of the 
host galaxy may be present, but it is also possible that systematic uncertainties in the SFR 
determinations are responsible. A similar study with a larger sample is needed to establish 
whether or not weak trends exist. 

Our results are consistent with previous findings that the luminosity of the brightest 
cluster in a galaxy is approximately proportional both to the number of clusters in the 
galaxy (e.g., Whitmore 2003) and to its star formation rate (Larsen 2002; Bastian 2008). 
However, inferences based only on the brightest cluster in a galaxy require extra (though 
plausible) assumptions about all the other, fainter clusters: their mass function, mass-to- 
light ratios, interstellar extinctions, and so forth. The method presented here is much more 
robust because it is based on the mass functions determined from hundreds or thousands of 
clusters over a wide range of masses and specific, narrow ranges of age (r < 10 Myr and 
100-400 Myr). 

Our results are also consistent with our previous findings that the mass and age distri¬ 
butions of clusters are similar in different galaxies. We approximate the joint distribution by 
a double power law: g{M, r) oc M^t'^ with (3^—2 and —1.0 ^ 7 ^ —0.5 (Fall & Chandar 
2012, Chandar et ah 2014). We refer to this as a quasi-universal distribution because the 
observed exponents differ relatively little among galaxies (typically, A/3 rs-/ Ay ~ 0.2)@ Our 
results for the CMF/SFR relation imply that the populations of clusters on galaxy scales 
are nearly the same apart from an overall normalization factor. The results presented here 
indicate that this normalization factor is the star formation rate. 

Our claim that populations of clusters are similar in different galaxies might seem to be 
contradicted by the fact that massive clusters are found in some galaxies but not in others. 
For example, the most massive clusters in the Antennae have Mmax ~ 10 ^ Mq whereas 
those in the SMC have Mmax ~ 10'^ M© (see Figure 1). However, this is largely, if not 
entirely, due to the different sizes of these samples or populations, roughly 10^ times larger 
in the Antennae than in the SMC. Indeed, it is easy to show that, for the quasi-universal 
model, the maximum cluster mass scales with the population size and hence with the star 


^ While the mass and age distributions of star clusters are remarkably similar on galaxy scales, in small 
regions there can be large variations in the formation and disruption rates, and hence in the observed mass 
and age distributions of the clusters (see e.g., Bastian et al. 2012; Silva-Villa et al. 2014; Chandar et al. 
2014 for an example of this effect in M83). As these regions are combined into larger ones, the variations 
will average out and the differences diminish. In particular, the exponent 7 of the age distribution may vary 
significantly with galactocentric distance for the physical reasons discussed by Fall & Chandar (2012). 
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formation rate as Mmax oc SFR where f5 again is the exponent of the mass function; 

thus Mmax oc SFR for j3 = —2. 

The quasi-universal model implies that the clusters in different galaxies are formed 
and disrupted in similar (although not exactly the same) ways. This is a remarkable fact, 
probably the consequence of a great deal of similarity in the structure of the interstellar 
media among galaxies on the scale of molecular clumps and clouds, in particular in their 
mass functions, which play roles in both the formation and disruption of clusters (see Fall, 
Krumholz, & Matzner 2010 and Fall & Chandar 2012 for further discussion of these topics). 
Until recently, it has not been possible to observe the ISM of other galaxies with sufficient 
angular resolution to test this conjecture, but with the advent of ALMA, this may now be 
within reach. 
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Fig. 1.— Mass functions of star clusters younger than 10 Myr in 7 star-forming galaxies. 
The symbols are as follows: LMC: pink; SMC: blue; NGC 4214: red; NGC 4449: orange; 
M83: green; M51: gray; and Antennae: purple. The selection of the clusters is described in 
the text. The left panel shows the unnormalized mass functions. The right panel shows the 
mass functions normalized by the current star formation rate (SFR). The fractional coverage 
of each galaxy is also accounted for. 
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Fig. 2.— Same as Figure 1, but now for clusters in the age range 100-400 Myr. 



Fig. 3.— Illustration of our method for determining the scatter in the CMF/SFR relation 
for young clusters (r < 10 Myr, left panel) and intermediate-age clusters (r = 100-400 Myr, 
right panel). The standard deviations of logR are given in the upper right of each panel 
(denoted by ua). 
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Fig. 4.— Logarithmic residuals in the CMF/SFR relations for young (r < 10 Myr, Filed 
circles) and intermediate-age (100-400 Myr, open circles) clusters plotted against a number 
of galaxy properties. No statistically signiFcant correlations are found between the residuals 
and the galaxy properties considered here for the young r < 10 Myr clusters. There are 
weak trends with galactic SFR, luminosity, and B — V color for the intermediate-age (100- 
400 Myr) clusters, but only the one with B — V is statistically signiFcant (at the 95% level). 
The top-left panel shows that our results do not agree with the two correlations presented 
in Kruijssen 2012, which were based on a compilation of observational results (dashed line) 
and theoretical calculations (dotted line). See text for more details. 




Table 1. Properties of Sample Galaxies 


Galaxy 

Name 

Galaxy 

Type 

Distance 

(Mpc) 

/S' 

r < 10 Myr 

Ha ± 25/im 
SFR2 

(A/e/yr) 

Published 

Ha SFR 
(Aifo/yr) 

Published 

Other SFR 
(AAe/yr) 

SFR 

References 

LMG 

Irregular 

0.05 

-1.60 ±0.10 

0.25 

0.25 

0.39 

1,2 

SMG 

Irregular 

0.06 

-1.69 ±0.15 

0.06 

0.04 

0.05 

3,4 

NGG 4214 

Irregular 

3.1 

-2.28 ±0.17 

0.11 

0.16 

0.22 

5, 5 

NGG 4449 

Irregular 

3.8 

-1.68 ±0.21 

0.35 

0.66 

0.89 

5, 5 

M83 

Spiral 

4.5 

-1.95 ±0.12 

2.65 

3.35 

5.01 

5, 5 

M51 

Spiral 

8.2 

-2.06 ±0.05 

3.20 

4.48 

7.58 

5, 5 

Antennae 

Merging 

22.2 

-2.18 ±0.04 

20.20 

25.48 


6 


^Power-law index for the mass function of star clusters, determined from the least-squares hts to 
log((iiV/(ilogM) = (/5 — 1) logM -|- const 

^Star formation rates determined in this work. See Section 3.1 for a description 


Note. — References: 1. James et al. 2008; 2. Harris & Zaritsky 2009; 3. Bolatto et al. 2011; 4. Harris 
& Zaritsky 2004; 5. Lee et al. 2009; 6. Zhang et al. 2001 
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Table 2. Uncertainties in SFR Estimates 


Source 

of Uncertainty 

Uncertainty 
log SFR 

Cluster Mass Function 

0.1 (r < 10 Myr) 

0.04 (r = 100-400 Myr) 

Flux Measurements 

~0.04-0.10 

Extinction Correction 

^0.1 

Fractional Coverage 

P30.04 

Lyc Escape Fraction 

unknown 

Metallicity 

systematic trend 
^0.4 at 1/10 Z© 


Table 3. Spearman Correlation Coefficients for Residuals Shown in Figure 4 


Galaxy 

Property 

Age Range 

rs 

r < 10 Myr 
P(> rs) 

Age Range r 

rs 

= 100-400 Myr 
P(> rs) 

log SFR (Me/yr) 

0.39 

0.39 

0.71 

0.07 

Ms (mag) 

-0.11 

0.82 

-0.66 

0.10 

B — V (mag) 

-0.29 

0.53 

0.78 

0.04 

inclination (deg) 

0.39 

0.38 

-0.50 

0.25 

Max gas velocity (km/s) 

0.25 

0.59 

0.61 

0.15 

CMF index {3 

0.07 

0.88 

-0.04 

0.94 


Note. — The nonparametric Spearman correlation coefficient, rs and the corre¬ 
sponding probabilities p(> rs) are compiled. A perfect correlation would have a 
value of ±1, and the sign indicates the direction of the association between the listed 
galaxy properties and the residuals in log A. Values for p(> rs) should be less than 
0.05 to have greater than 95% conhdence that there is a real correlation. The only 
correlation that satishes this condition is the one with B — V for intermediate-age 
clusters. 












